Efficient Mixing at low Reynolds numbers using polymer additives 
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Mixing in fluids is a rapidly developing fleld of fluid mechanics Jll-Q, being an 
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important industrial and environmental problem. The mixing of liquids at low 
Reynolds numbers is usually quite weak in simple flows, and it requires special 
devices to be eflicient. Recently, the problem of mixing was solved analytically 
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Ph \ for a simple case of random flow, known as the Batchelor regime @-|^. Here we 



demonstrate experimentally that very viscous liquids at low Reynolds number. 
Re. Here we show that very viscous liquids containing a small amount of high 



^ ' molecular weight polymers can be mixed quite efliciently at very low Reynolds 
^ ' numbers, for a simple flow in a curved channel. A polymer concentration of only 



0.001% suflices. The presence of the polymers leads to an elastic instability [^] 



and to irregular flow [ |10|| , with velocity spectra corresponding to the Batchelor 



C ■ regime [^j-]^. Our detailed observations of the mixing in this regime enable us to 
> 

' conflrm sevearl important theoretical predictions: the probability distributions 

■ of the concentration exhibit exponential tails [^,^, moments of the distribution 

decay exponentially along the flow [Q, and the spatial correlation function of 

concentration decays logarithmically. 

Solutions of flexible high molecular weight polymers are viscoelastic fluids |]rT|. There 
are elastic stresses that appear in the polymer solutions in a flow, and that grow non-linearly 



with the flow rate [ITJ]. This can lead to many special flow effects, including purely elastic 



transitions 0,|T2|,[T3| that qualitatively change character of the flow at vanishingly small Re. 



As a result of such transitions secondary vortical flows appear in different systems, where 



the primary motion is a curvilinear shear flow. Onset of those secondary flows depends on 
the Weissenberg number, Wi = A7, where A is the polymer relaxation time and 7 is the 
shear rate. It plays a role analogous to that of Re in competition between non-linearity 
and dissipation. As it has been reported recently, an elastic flow transition can result in a 
special kind of turbulent motion, elastic turbulence which arises at arbitrary small Re. 
In our experiments we flnd that irregular flow excited by the elastic stresses in the polymer 
solutions can lead to quite efficient mixing at very low Re. 

Our mixing experiments were carried out in a curvilinear channel. Fig. la, at a room 
temperature, 22.5±0.5 °C. The total rate of the fluid discharge, Q, was always kept constant, 
so that the average time of mixing was proportional to the position, N, along the channel. 
Fig. la. We used a solution of 65% saccharose and 1% NaCl in water, viscosity rjg = 0.153 
Pa s, density p = 1.32 g/cm^, as a solvent for the polymer. We added polyacrylamide, 
Myj = 18, 000, 000, at a concentration of 80 p. p.m. by weight. One of the solutions entering 
the channel. Fig. la, also contained Cq = 2 p. p.m. of a fluorescent dye. The solution viscosity 
was f] = 0.198 Pa-s at a shear rate 7 = 4 s~^. The relaxation time. A, estimated from the 
phase shift between the stress and the shear rate in oscillatory tests, was 1.4 s. An estimate 
for the diffusion coefficient of the dye is given by that for the saccharose molecules, which is 
about D = 8.5 ■ 10~^ cm^/s. The characteristic shear rate and the Weissenberg number in 
the flow are estimated as 7 = {2Q / (f ) / {d/2) = AQ/d? and Wi = \{AQ/d?), respectively. 

The Reynolds number. Re = 2Qp/{dri), was always quite small. It only reached 0.6 for 
highest Q that we explored. Therefore, flow of the pure solvent remained quite laminar 
and no mixing occurred. Fig. lb. The boundary separating the liquid with and without 
the dye was smooth and parallel to the direction of the flow, and it only became smeared 
due to molecular diffusion as the liquids advanced downstream. Behaviour of the polymer 
solution was qualitatively different from that of the solvent. The flow only was laminar and 
stationary up to a value of Q corresponding to Wic = 3.2 (and Re = 0.06), at which an 
elastic instability occurred. This instability lead to irregular flow and mixing of the liquids, 
Fig.lc. The experiments were carried out at Q about twice above the flow instability onset. 



Wi = 6.7, at which homogeneity of the mixture at the exit of the channel was the highest. 

Power spectra of fluctuations of velocity in the center of the channel at Wi = 6.7 are 
shown in Fig.2. The spectra of both longitudinal and transversal velocity components do 
not exhibit any distinct peaks and have broad regions of a power law decay, that is typical 
for turbulent flow. They resemble the power spectra found at similar Wi in flow of the same 
polymer solution between two parallel plates in the regime of elastic turbulence [l^. So, we 
believe the origin of the irregular flow is the same here as in Ref. [|10| . 

According to the Taylor hypothesis fluctuations of the flow velocity in time are mainly 
due to fluctuations in space advected by the mean flow So, the spectra in Fig.2 imply 
that the power of the velocity fluctuations scales with the fc-number in space as P ~ k~^'^. 
Fluctuations of the velocity gradients should scale then as k~^'^, so that the flow becomes 
increasingly homogeneous at small scales, and mixing is mainly due to the largest eddies 
having the size of the whole flow system 0. Such flow conditions correspond to the Batchelor 
regime of mixing 1^0. It is one of the two simple cases of random flow, where the problem 
of mixing has been solved analytically Investigation of mixing in the irregular flow 

in the channel gave us an opportunity to study mixing in the Batchelor regime in detail and 
to test a few key the theoretical predictions for the first time. 

Mixing of the polymer solution was a random process. So, it was appropriate to charac- 
terize it statistically, by probability distribution function (PDF) to find different concentra- 
tions, c, of the dye in a point, and by values of the moments, Mj, of the distribution. An 
i-th moment is defined as an average, < |c — c|* > /<?, where c is the average concentration 
of the dye, which in our case is equal to Cq/2. Small values of Mj mean high homogeneity 
and good mixing of the liquids. Mj are all unity at the channel entrance and they become 
zero for perfectly mixed liquids. 

PDF of the concentrations at = 8 and A^ = 24 (with Mi = 0.25 and 0.72, respectively) 
are shown in Fig. 3a-b. Dependencies of Mi and M2 on the position, A^, along the channel 
are presented in Fig.4. Advanced stages of mixing corresponding to A^ > 30 were studied, 
when the liquids were premixed before they entered the channel, Fig.4. Representative PDF 



in the flow with prehminary mixing taken at the positions, where Mi = 0.082 and 0.030, 
are shown in Fig.3c-d, respectively. 

Representative spatial autocorrelation functions for the dye concentration are shown in 
Fig. 5. At the beginning the concentration distribution is strongly influenced by the initial 
conditions. There are extended uniform regions with maximal and zero dye concentration, 
so, the concentration remains highly correlated over rather large distances. Fig. 5a, and PDF 
has maxima near cq and zero, Fig.Sa. As the liquid advances downstream, it becomes in- 
creasingly homogeneous. Mixing patterns exhibit a lot of flne scale structures of different 
brightnesses, the correlation distances become shorter. Fig. 5b, and PDF of the dye concen- 
tration become narrower, Fig.Sb. PDF in Fig.3b has a single peak at c and long tails that 
decay exponentially and touch the limits of the concentration, zero and Cq. 

Further downstream mixing patterns have characteristic features at similar spatial scales, 
Fig.5c-e, but are much more faded. The PDF in Fig. 3c is much narrower than in Fig. 3b and 
has quite clear exponential tails, in agreement with the theoretical predictions that 
implies signiflcant intermittency in mixing 0. The distribution is well conflned in a region 
far from the limits of zero and Cq. So, the dependence on the initial conditions should be 
quite minor by that point. At the last point (A^ = 54, Mi = 0.030) the non-homogeneity in 
the concentration is hardly seen, and the PDF, Fig. 3d, is very narrow. One can see that at 
A^ of about 30 the autocorrelation function reaches some asymptotic form, which does not 
change as the liquid moves further downstream. This may correspond to a situation, when 
the small scale structures are smeared by the molecular diffusion at the same rate as they are 
created by the flow. The autocorrelation functions, Fig.5c-e, decay logarithmically at the 
distances above the diffusion length, that is in agreement with the theoretical predictions for 
the Batchelor regime [Q-Q. Spatial structure of concentration fluctuations in the Batchelor 
regime has been experimentally studied before in different types of flows 0,0,^ . However, 



agreement between experiment and theory remained rather controversial (see Ref.2, 14 and 
15 for discussion). 

The dependencies of Mi and M2 on A^ in Fig.4 have inflections at A^ between 20 and 30. 



From the above discussion we can suggest that they are due to transition to an asymptotic 
regime, where dependence on the initial conditions is lost. The correlation functions at 
different positions become identical, and PDF of concentrations can be superimposed rather 
well for Ac of about 3Mic around c by appropriate stretching of the axes. There is no self- 
similarity in the shapes of PDF at larger Ac/(Mic), however, and the shape permanently 
changes with A^. One can learn from Fig. 4 that both Mi and M2 decay exponentially above 
= 30, the rate of the decay being two times higher for M2 than for Mi. The higher 
order moments were found to decay exponentially, Mj ~ exp{—'~fiN), as well, that is quite in 
agreement with the theoretical predictions P]. We note here that M2, which is the variance 
of fluctuations of the dye concentration, is often considered as an analogue of energy of 
turbulent motion 0. So its logarithmic derivative, ^^^^ = —72/7.8 s^^ (see Fig. 4), is 
the analogue of the energy dissipation rate. The theory also predicts linear growth of the 
coefficients, 7^, with i at small i and saturation of the growth, when i becomes large 0. This 
form of dependence of 7^ on i is quite confirmed by our experimental results, inset in Fig.4. 
The deviation of the dependence of 7^ on i from a straight line at large i is a quantitative 
evidence for the lack of self-similarity in the mixing. One can see that all our experimental 
results agree very well with all the theoretical predictions for the Batchelor regime of mixing, 
we tested ||-||. 

As it is suggested by the decay of Mi in Fig.4, the polymer solutions are mixed at a 
characteristic distance of AA^ ~ 15. It corresponds to the total path of about lAOd and the 
average fiow time of about 120 seconds, which is three orders of magnitude smaller than the 
diffusion time, d?/D. Using a more concentrated sugar syrup as a solvent, we prepared a 
polymer solution with D about 30% smaller and with viscosity and relaxation time about 
two times larger than those of the original solution. We studied distribution of c at A^ = 29 in 
a flow of this solution at the same Wi of 6.7 as before. Here the characteristic flow velocities 
were twice lower. Re was about four times lower, and the ratio between the flow time and the 
diffusion time, (P/D, was about 1.5 higher than in the original flow. The measured values 
of Ml and M2 were the same as in Fig.4. So, the inertial forces and the molecular diffusion 
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did not have any apparent influence on the mixing efficiency. Dependence of the efficiency 
of mixing at the optimal flow conditions on concentration of the polymers was surprisingly 
weak (although Wic grew fast, when the polymer concentration was decreasing). So, for a 
solution with the polymer concentration of 10 p. p.m. (77/775 = 1.03), Mi of as low as 0.22 
was reached at A?" = 29 (and at Re — 0.065). The mixing was observed down to the polymer 
concentration of 7 p. p.m. Thus, very viscous liquids can be efficiently mixed in curvilinear 
channels at very low flow rates with the aid of polymer additives at very low concentrations. 
This method of mixing, we believe, can find some industrial and laboratorial applications. 
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FIGURES 

FIG. 1. Experimental set up and two snapshots of the flow, a) Experimental set-up. A 
channel of a depth d = 3 mm is machined in a transparent bar of perspex and sealed from above 
by a transparent window. The channel is square in the cross-section, and consists of a sequence 
of 60 smoothly connected half-rings with the inner and outer radii -Ri = 3 mm and i?2 = 6 mm, 
respectively. The flow is always observed near the middle of a half-ring on the right-hand side of 
the bar, when viewed streams. The half-rings on the right-hand side are numbered starting from 
the channel entrance, and the number, N, of the half-ring is a natural linear coordinate along 
the channel. Two liquids are fed into the channel by two syringe pumps at equal discharge rates 
through two separate tubes. The two working liquids are always identical with the only difference of 
a small amount of a fluorescent dye (fluorescein) added to one of them. The channel is illuminated 
from a side by an Argon-Ion laser beam converted by two cylindrical lenses to a broad sheet of 
light, with a thickness of about 40 /xm in the region of observation. The fluorescent light emitted 
by the liquid in the perpendicular direction is projected onto a CCD camera and digitized by a 
8-bit, 512x512 frame grabber. Concentration of the dye is evaluated from the intensity of the light, 
which was found to be proportional to the concentration, b-c Snapshots of the flow at A'^ = 29. 
Bright regions correspond to high concentration of the fluorescent dye. Curved boundaries of the 
channel are seen on the left and on the right, b) Pure solvent at Re = 0.16. c) Polymer solution 
at the same flow rate, corresponding to Wi = 6.7. 

FIG. 2. Power, P, of fluctuations of velocity in the middle of the channel at A" = 12 as a 
function of frequency, /. The flow velocity was measured by a laser Doppler anemometer, when 
the region of the laser beam crossing was made very small, 15 x 15 x 40 //m, to decrease the gradient 
noise. The spectra of the velocity components for the polymer solution along the mean flow, curve 
a, across the mean flow, curve b, and for the pure solvent across the mean flow, curve c, are shown. 
The mean velocity for the polymer solution was V = 6.6 mm/s. The RMS of fluctuations, Vrms, 
was 0.09y and 0.04y in the longitudinal and transversal directions, respectively. 
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FIG. 3. Plots of PDF of concentration of the fluorescent dye at different positions. Brightness 
profile was taken 12.5 times per second along a single line across the channel near the middle of 
a half-ring (a horizontal line in the middle in Fig.lb-c). Each plot represents statistics over about 
10^ points, corresponding to about 25,000 different moments of time, and a total liquid discharge 
of 2 • 10^ (f. The regions near the walls of the channel with the widths of O.ld were excluded from 
the statistics, a — h- PDF at = 8 and = 24; c — d - with preliminary mixing, at the positions 
corresponding to AT = 39 and N = 54, respectively. 

FIG. 4. Dependence of Mi and M2, circles and squares, respectively, on the position. A'', along 
the channel. The average flow time is connected to N as t = N ■ 7.8 s. In order to observe stages 
of mixing corresponding to N > 30, we carried out a series of experiments, where the liquids were 
premixed before they entered the channel. A shorter channel of the same shape was designed 
for this purpose and put before the entrance to the original one. As a result of this preliminary 
mixing, PDF of the dye concentrations at N = 2 was almost identical to PDF at N = 27 without 
the premixing. Filled symbols represent Mi and M2 without preliminary mixing. Empty symbols 
are for the flow with preliminary mixing, when a value of 25 is added to the actual position, in 
order to match the entrance conditions. . One can see that these curves are indeed continuations 
of the dependencies obtained for Mi and M2 in the channel without the premixing. 

FIG. 5. Correlation coefficients for the concentration as functions of the distance. Ax, across 
the channel (semilogarithmic coordinates). The curves a — d correspond to the same conditions 
as PDF in Fig.5a-d, and they were calculated using the same data arrays. The curve e is for 
A'^ = 29. The correlation functions coincide for A^ > 29. At small Ax they have parabolic scaling 
with a characteristic length xq = O.OlTd 50 /xm. It is, probably, defined by the thickness of the 
illuminating light sheet, about 40 /xm, and by the molecular diffusion scale, {D /{Vrms/d))^/'^ ~ 25 
jjLm.. At larger Ax the scaling is logarithmic. 
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